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Summary 



The emergence of a new generation of instrumentation in astrophysics, which provide spatially-resolved 
spectra over a large 2-dimensional (2D) field of view, offers the opportunity to perform emission-line surveys 
based on samples of hundreds of spectra in a 2D context, enabling us to test, confirm, and extend the previous 
body of results from small-sample studies based on typical long-slit spectroscopy, while at the same time 
opening up a new frontier of studying the 2D structure of physical and chemical properties of the disks of 
nearby spiral galaxies. The project developed in this dissertation represents the first endeavour to obtain 
full 2D coverage of the disks of a sample of spiral galaxies in the nearby universe, by the application of 
the Integral Field Spectroscopy (IFS) technique. The semi-continuous coverage spectra provided by this 
spectral imaging technique allows to study the small and intermediate linear scale variation in line emission 
and the gas chemistry in the whole surface of a spiral galaxy. 

The PPAKIFS Nearby Galaxies Survey: PINGS, was a carefully devised observational project, designed 
to construct 2D spectroscopic mosaics of 17 nearby galaxies in the optical wavelength range. The sample 
includes different galaxy types, including normal, lopsided, interacting and ban^ed spirals with a good range 
of galactic properties and star forming environments, with multi-wavelength public data. The spectroscopic 
data set comprises more than 50 000 individual spectra, covering an observed area of nearly 100 arcmin^, 
an observed surface without precedents by an IFS study. All sources of eiTors and uncertainties during the 
reduction process of the IFS observations are assessed very carefully. This methodology contributed not only 
to improve the standard reduction pipeline procedure for the particularly used instrument, improvements that 
can be applied to any similar integral-field observation and/or data reduction, but to defining a self-consistent 
methodology in terms of observation, data reduction and analysis for the kind of IFS surveys presented in 
this dissertation, as well as providing a whole new set of IFS visualization and analysis software made 
available for the public domain. 

The scientific analysis of this dissertation comprises the study of the integrated properties of the ionized 
gas of the whole PINGS sample, and a detailed 2D study of the physical and chemical abundance distribu- 
tion derived from the emission line spectra of four selected galaxies of the sample. Spatially-resolved maps 
of the emission line intensities and physical properties are derived for each the selected galaxies. Different 
methodologies are explored in order to study the spatially-resolved spectroscopic properties of the galaxies. 
Abundance analysis are performed based on a variety of diagnostic techniques using reddening corrected 
spectra. From this analysis, evidence is found to support that the measurements of emission lines of a "clas- 
sical" H II region are not only aperture, but spatial dependent, and therefore, the derived physical parameters 
and metallicity content may significantly depend on the morphology of the region, on the slit/fibre position, 
on the extraction aperture and on the signal-to-noise of the observed spectrum. On the other hand, the results 
presented in this dissertation indicate the existence of non-linear multi-modal abundance gradients in normal 
spiral galaxies, consistent with a flattening in the innermost and outermost parts of the galactic discs, with 
important implications in terms of the chemical evolution of galaxies. 

The powerful capabilities of wide-field 2D spectroscopic studies are proven. The chemical composition 
of the whole surface of a spiral galaxy is characterised for the first time as a function not only of radius, but of 
the intrinsic morphology of the galaxy, allowing a more realistic determination of their physical properties. 
The methodology, analysis and results of this dissertation will hopefully contribute in a significant way to 
understand the nature of the physical and chemical properties of the gas phase in spiral galaxies. 

ill 



Declaration 



I hereby declare that this thesis entitled Multi-dimensional analysis of the chemical and physical properties 
of spiral galaxies is not substantially the same as any that I have submitted for a degree or diploma or other 
qualification at any other University. I further state that no part of this dissertation has already been or is 
being concurrently submitted for any such degree, diploma or other qualification. This thesis is essentially 
the result of my own work, and includes nothing which is the outcome of work done in collaboration except 
where specifically indicated. Those parts of this dissertation which are undergoing review for publication, 
or included in conference proceedings are as follows: 

• Chapter 3 

The work presented in this chapter is essentially my own, it has been submitted for publication as: 
Rosales-Ortega, F. F., Kennicutt, R. C, Sanchez, S. F, Diaz, A. I., Pasquali, A., Johnson, B. D. and 
Hao, C. N., (2009) PINGS: the PPAK IF S Nearby Galaxies Survey, submitted to Monthly Notices of 
the Royal Astronomical Society, and it was benefited from collaboration with these authors. 

• Chapter 4 

The work presented in this chapter is essentially my own, a large part of it was submitted for publi- 
cation in the article mentioned above. The work related to NGC 628 has been partially done within 
a collaboration, and submitted for publication as: Sanchez, S. F, Rosales-Ortega, F. F, Kennicutt, 
R. C, Diaz, A. I., Pasquah, A., Johnson, B. D. and Hao, C. N., (2009) PPAK Wide-field Integral Field 
Spectroscopy of NGC 628: The largest spectroscopic survey on a single galaxy, submitted to Monthly 
Notices of the Royal Astronomical Society. The Gaussian-suppression and absolute fux calibration 
techniques were developed in collaboration with S. F. Sanchez. Their implementation and initial test- 
ing were carried out in parallel by both of us, though the final implementation for this dissertation are 
my own work. 

• Chapter 5 

The work presented in this chapter is essentially my own. Some parts has been partially done within 
a collaboration, and submitted for publication in the articles mentioned above. 

• Chapter 6 

The work presented in this chapter is essentially my own, but benefited from discussions and advice 
from S. F. Sanchez and A. I. Diaz. 

• Some figures of chapter 3 and chapter 6 have been included in The Promise of Multiwavelength 
and IFU observations, Kennicutt, R. C, Hao, C. N., Johnson, B. D., Rosales-Ortega, F. F, Diaz, 
A. I., Pasquali, A. and Sanchez, S. F, (2009) Proceedings of the lAU Symposium 262, G. Bruzual, 
S. Chariot, eds. 

This thesis is less than 60 000 words in length. 

Fernando Fabian Rosales Ortega 
Cambridge, November 5, 2009 



V 



Thesis content 



The whole thesis is not included in astro-ph due to file size limitations. 

The full contents can be found at: 
http : / / www . dspace . cam . ac . uk/handle/ 1810/224843 



vii 



1 

Introduction 




he existence and distribution of the chemical elements and their isotopes in the universe is a 
consequence of very complex processes that have taken place in the past since the Big Bang and 
subsequently in stars and in the interstellar medium (ISM) of the present day galaxies, where 
they are still ongoing. These processes have been studied theoretically, experimentally and ob- 



servationally. Different theories of cosmology, stellar evolution and interstellar processes have been consid- 
ered, laboratory investigations of nuclear and particle physics, studies of elemental and isotopic abundances 
in the Earth and meteorites have also been involved, as well as astronomical observations of the physical 
nature and chemical composition of stars, galaxies and the interstellar medium. 

From the observational point of view, the study of chemical abundances in galaxies, like many other 
areas of astrophysics, has undergone a remarkable acceleration in the flow of data over the last few years. 
We have witnessed wholesale abundances determinations in tens of thousands of galaxies from large scale 
surveys such as the Two Degree Field Galaxy Redshift Survey (2dFGRS, CoUess et al., 2001) and the Sloan 
Digital Sky Survey (SDSS, York et al., 2000), measurements of abundances in individual stars of Local 
Group galaxies beyond the immediate vicinity of the Milky Way, and the determination of the chemical 
composition of some of the first stars to form in the Galactic halo. Chemical abundances studies are also 
increasingly being extended to high redshift, charting the progress of stellar nucleosynthesis over most of 
the age of the universe. The primary motivation common to all of these observational efforts is to use the 
chemical information as one of the means at our disposal to link the properties of high redshift galaxies with 
those we see around us today, and thereby understand the physical processes at play in the formation and 
evolution of galaxies. 

The galactic chemical evolution is dictated by a complex array of parameters, including the local initial 
composition, star formation history (SFH), gas infall and outflows, radial transport and mixing of gas within 
disks, stellar- yields, and the initial mass function (IMF). Although is difficult to disentangle the effects of 
the various contributors, measurements of current elemental abundances constrain the possible evolutionary 
histories of the existing stars and galaxies. Important constraints on theories of galactic chemical evolution 
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and on the star formation histories of galaxies can be derived from the accurate determination of chemical 
abundances either in individual star-forming regions distributed across galaxies or through the comparison of 
abundances between galaxies. Nebular emission lines from individual H II regions have been, historically, 
the main tool at our disposal for the direct measurement of the gas-phase abundance at discrete spatial 
positions in low redshift galaxies. 

However, in order to obtain a deeper insight of the mechanisms that rule the chemical evolution of galax- 
ies, we require a significantly the number of H II regions sampled in any given galaxy. In this dissertation, 
I present a new observational technique conceived to tackle the problem of the 2-dimensional coverage of 
the whole surface of a galaxy. The advent of new spectroscopic techniques provides powerful tools for 
studying the small and intermediate scale-size variation in line emission and stellar continuum in nearby 
well-resolved galaxies. In this work, I address the problems and challenges that imply the determination 
of the chemical composition in galaxies in a 2D context and the subsequent derivation of their physical 
properties. 

I will begin by presenting in this chapter a literature review on the determination of chemical abundances 
in galaxies. As an introduction to this topic, the physics of gaseous nebulae is discussed in § 1.1, together 
with a discussion of extra-galactic Hll regions in § 1.2. The different methods of abundance determinations 
are presented in § 1.3. Physical properties derived from the determination of chemical abundances are dis- 
cussed in § 1 .4. These latter sections ai^e partially based on the paper reviews and books about the physics 
and chemistry of the interstellar medium and H II extragalactic regions by Dinerstein (1990), Perez-Montero 
& Diaz (2005), Tielens (2005), and Osterbrock & Ferland (2006). This discussion leads to the presenta- 
tion of new techniques and methods for the determination of chemical abundances in nearby galaxies as 
described in § 1.5 

1.1 The Physics of Gaseous Nebulae 

Gaseous Nebulae are observed as bright extended objects in the sky, some are easily observed on direct 
images but many others are intrinsically less luminous or are affected by interstellar extinction on ordinary 
images, but can be resolved on long exposures with special filters and techniques so that the background and 
foreground stellar and sky radiations are suppressed. The surface brightness of a nebula is independent of its 
distance, but more distant nebulae have (on average) smaller angular size and greater interstellar extinction. 

Gaseous nebula have emission-line spectra. Nebulae emit electromagnetic radiation over a broad spectral 
range, although only a few wavelengths pass easily through the Earth's atmosphere. Visible light and some 
infrared and radio radiation can be studied from the ground, but most other wavelengths can only be covered 
from high-latitude aircrafts or space telescopes. The source of energy that enables normal emission nebulae 
to radiate is ultraviolet radiation from stars within or near the nebula. There should be one or more stars 
with effective surface temperature T^, > 3 x 10"^ K, the ultraviolet photons of these stai^s transfer energy 
to the nebula by photoionization. In all nebulae, hydrogen (H) is by far the most abundant element, and 
photoionization of H is thus the main energy-input mechanism. Photons with energy greater than 13.6 
eV (the ionization potential of H), are absorbed in the process, and the excess energy of each absorbed 
photon over the ionization potential appears as kinetic energy of a newly liberated photoelectron. Collisions 
between electrons and between electrons and ions, distribute this energy and maintain a Maxwellian velocity 
distribution with temperature T in the range 5,000 < T < 20,000 K in typical nebulae. Collisions between 
thermal electrons and ions excite the low-lying energy levels of the ions. Downward radiation transitions 
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from these excited levels have very small transition probabilities, but at the low densities (Ue < 10^ cm^^) of 
typical nebulae, collisional de-excitation is even less probable, so almost every excitation leads to emission 
of a photon, and the nebula thus emits a forbidden-line spectrum that is quite difficult to reproduce under 
terrestrial laboratory conditions. 

Thermal electrons ai^e recaptured by the ions, and the degree of ionization at each point in the nebula is 
fixed by the equilibrium between photoionization and recapture. In the recombination process, recaptures 
occur to excited levels, and the excited atoms thus formed then decay to lower and lower levels by radiative 
transitions, eventually ending in the ground level. In this process, line photons are emitted and this is the 
origin of the Hi B aimer and Paschen line spectra observed in all gaseous nebulae. The recombination 
of H"*" gives rise to excited atoms of H'' and thus leads to the emission of the H I spectrum. Likewise, He"'" 
recombines and emits the He I spectrum, and in the most highly ionized regions, He++ recombines and emits 
the He II spectrum. Recombination lines of trace elements are also emitted; however, the main excitation 
process responsible for the observed strengths of such lines with the same spin or multiplicity as the ground 
term is resonance fluorescence by photons, which is much less effective for H and He lines because the 
resonance lines of these more abundant elements have greater optical depth. Nevertheless, line emission of 
these rare elements plays a significant role in the physics of the nebula, and permits the determination of the 
chemical composition inside the nebula. 

The spectra of gaseous nebulae are dominated by coUisionally excited forbidden lines of ions of com- 
mon elements, such as [OlII] AA4959, 5007 (the famous green nebular lines); [Nil] AA6548, 6584 and 
[Sll] AA9069,9523 in the red; and [Oil] AA3727, 3729 in the ultraviolet (which normally appears as a 
blended A3727 line on low-dispersion spectrograms). In addition, the permitted lines of hydrogen. Ha 
A6563 in the red, HjS A4861 in the blue, Hy A4340 in the violet and so on, are characteristic features of 
every nebular spectrum, as is He I A5876, which is considerably weaker, while Hell A4686 occurs only 
in higher-ionization nebulae. Long-exposure spectrophotometric observations extending to faint intensities 
show progressively weaker forbidden lines, as well as faint permitted lines of common elements such as 
C II, C III, C IV and so on. Nebular emission-line spectra extend into other spectral ranges, in the infrared for 
example, the [Nell] A 12.8 /im and [OlII] A88.4 /xm are among the strongest lines measured, into the ultra- 
violet Mgll AA2796,2803, Cm] A A 1907, 1909, Civ A A 1548, 1551 and Lya A 1216 are also observed. Is 
often necessary to obtain spectra outside the traditional visible/near spectral bands to get an accurate picture 
of the system in question. 

Gaseous nebulae have weak continuous spectra, consisting of atomic and reflection components. The 
atomic continuum is emitted chiefly by free-bound transitions, mainly in the Paschen continuum of H I at 
A > 3646 A, and the Bahner continuum at9I2A <A< 3646 A. In addition to the bright-line and continuous 
spectra emitted by atomic processes, many nebulae have reflection continua arising from starlight scattered 
by dust. The amount of dust varies from nebula to nebula, and the strength of this continuum fluctuates 
correspondingly. In the infrared for example, the nebular continuum is largely thermal radiation emitted by 
dust particles heated to a temperature of order 100 K by radiation derived originally from the central star. 

Gaseous nebula may be classified into two main types: H II regions and planetary nebulae (PNe). Though 
the physical processes in both types are quite similar, the two groups differ greatly in origin, mass, evolution 
and age. The objects of study for the chemical composition of galaxies are extragalactic H II regions. These 
diffuse nebulae are regions of interstellar- gas in which the exciting stars are O- or early B-type stars, i.e. 
young stai^s which use up their nuclear energy quickly. These hot, luminous stars undoubtedly formed fairly 
recent from interstellar matter that would be otherwise be part of the same nebula. The effective temperature 
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of the stars aie in the range 3 x 10^ < < 5 x 10^ K; throughout the nebula, H is ionized, He is singly 
ionized and other elements are mostly singly or doubly ionized. Typical densities in the ionized part of 
the nebula are of the order 10 or 10^ cm^^, ranging to high as 10^ cm^^. In many nebulae, dense neutral 
condensations are scattered through the ionized volume. Internal motions occur in the gas with velocities 
of order 10 kms^\ approximately the isothermal sound speed. Bright rims, knots, condensations, and so 
on, are apparent to the limit of resolution. The hot, ionized gas tends to expand into the cooler suiTounding 
neutral gas, thus decreasing the density within the nebula and increasing the ionized volume. The outer edge 
of the nebula is surrounded by ionization fronts running out into the neutral gas. This original two-phase 
model of the interstellar medium (the Hll /Hi region dichotomy) was introduced by Stromgren (1939). He 
showed that photoionized gas neai^ hot stars is segregated into physically distinct volumes, sepai^ated from 
their neutral environment by sharp boundaries. 

1.1.1 ExtragalacticH II regions 

The spectra of H II regions are strong in H I recombination lines, [N II] , [O II] and [O III] collisionally excited 
lines, but the strengths of N and O may differ greatly, being stronger in the nebulae with higher central- 
star temperatures. The brightest Hll regions can easily be seen on almost any large-scale image of nearby 
galaxies, and those taken in a narrow wavelength band in the red (including Ha and [N II] lines) are specially 
effective in showing faint and often heavily obscured extragalactic H II regions. The H II regions are strongly 
concentrated to the spiral arms and indeed are the best objects for tracing the chemical composition, structure 
and dynamics of the spiral arms in distant galaxies. They trace recent star formation and, thi^ough the analysis 
of their chemical composition, previous star formation activity. Typical masses of observed H II regions are 
of the order 10^ to 10"^ Mq, with the lower limit depending largely on the sensitivity of the observational 
method used. 

Hll regions are the only form of interstellar material which emits strongly in the optical spectral region; 
therefore, there is a much longer and richer history of observations and theory for them than for the other 
thermal phases of interstellar matter. Optical observations of H II regions provide fairly complete informa- 
tion about their elemental composition. From their spectra, abundances relative to hydrogen can be estimated 
for nearly all of the most common elements, particularly He, N, O, Ne, Ar, and S (note that oxygen alone 
constitutes nearly 50% by mass of the elements heavier than helium.) Furthermore, ionized nebulae are 
remarkably efficient machines for converting ultraviolet continuum energy from OB stars, originally diluted 
over wide bandpasses, into a few narrow, intense, optically-thin emission lines. The intrinsic emissivities 
of these lines aie easy to calculate in principle, although they are sensitive to the local thermodynamic state 
of the gas (electron density rig and temperature Tg). On the other hand, the thermal parameters can also be 
determined from the spectra, using diagnostic line-intensity ratios. In this way, Hll regions can be used to 
measure element abundances in the (present-day) gas of distant galaxies. 

The sample of extragalactic H II regions studied so far has metal abundances ranging from about 0.02 to 
several times solar. This is a useful complement to studies of our own Galaxy, which contains no severely 
metal-deficient H II regions (except for a handful of planetary nebulae formed by stars of the halo popu- 
lation). In contrast, for many Hll regions in the outskirts of late-type spirals and in some dwarf irregular 
galaxies, the process of metal enrichment by stellar nucleosynthesis is still in its eai^ly stages, providing a 
hint on the eaiiy chemical evolution of galaxies. These low-metallicity H II regions are also presumed to 
have experienced only a small degree of alteration in their hehum abundances due to stellar activity. There- 
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fore, their present He/H ratios should be nearly the same as the primordial value, providing valuable tests 
for cosmological theories. The various categories of extragalactic H II regions are essentially lists of their 
environments. These include: 

1. Disk H II regions in spiral and irregular galaxies. 

2. Gassy dwarf irregular- galaxies with spectra which are heavily dominated by H II regions. 

3. Nuclear and near-nuclear regions sometimes called "starburst" or "hotspot" Hll regions (e.g. Kenni- 
cutt et al., 1989). 

The first two categories have the best abundance data available in the literature. Regions in the third 
group tend to have relatively strong stellar continua and to be fairly metal-rich, which make it to difficult 
to obtain accurate measurements of the emission lines from which abundances are determined. On the 
other hand, members of the first two categories are universally regarded as members of the same family. 
Hll regions in nearby galaxies have been well-catalogued; atlases are available for the Large and Small 
Magellanic Clouds (LMC, SMC), and a large number of other galaxies (e.g. Hodge & Wright, 1967, 1977; 
Hodge & Kennicutt, 1983). The stai^-forming dwarf irregulars are usually found by spectroscopic surveys 
for emission-line galaxies (e.g. Kinman 1984 and more recently the SDSS data releases). 

The statistical properties of the H II region populations in spiral and irregular galaxies were addressed by 
Kennicutt (1988) and Kennicutt et al. (1989). They find that late-type galaxies have both intrinsically higher- 
luminosity Hll regions, and larger total numbers of Hll regions after normalization by galaxy size, than do 
eaily-type spirals. Within a galaxy, the differential luminosity function of the Hll regions is roughly a 
power-law, A'^ oc L^^^"^, although some low-luminosity irregulars have an exceptional supergiant complex, 
and Sa-Sb galaxies are deficient in luminous regions. While the positive correlation between the luminosity 
of the brightest H II region and that of the parent galaxy can be understood as chiefly a sample-size effect, 
the dependence on morphological type is a real and separate factor. Typical large galaxies contain hundreds 
of optically detectable H II regions. It is important to note that of all the regions detected and cataloged in 
Ha or H/3 , it is usually the nearest and the most luminous "giant" H II regions for which abundances are 
derived. 

Some of the best-studied regions are the 30 Doradus complex in the LMC, NGC 604 in M 33, and 
NGC 5461 and 5471 in M 101. Selection effects play an important role, necessarily poorer spatial resolution 
contributes to a tendency to identify larger regions in more distant galaxies. This effect is illustrated by Israel 
et al. (1975), who compare large-beam radio measurements with optical images of the same Hll regions in 
M 101; at better resolution these regions break up into groups or chains of smaller clumps. Likewise, H II re- 
gions in dwarf irregulars are also found to have complex structure when closely examined (e.g. Hodge et al., 
1989; Davidson et al., 1989). In more distant galaxies, we will always be looking at more heterogeneous 
volumes; for example, a typical aperture size (4") for spectrophotometric studies corresponds to 1 pc at 50 
kpc (the LMC) and 2 kpc at 100 Mpc. 

The morphology of many giant extragalactic H II regions can be characterized to first order as a "core- 
halo" structure, on the basis of both optical and radio-continuum data. The cores are composed of dense 
material, often in several distinct clumps, close to the ionizing stars. The diffuse, lower-density envelopes are 
presumably ionized by photons escaping from the inner regions and represent the radiation-bounded edges 
of the Sti"6mgren volume. Most giant extragalactic H II regions are believed to be essentially radiation- 
bounded (e.g. McCall et al, 1985). In addition, the denser regions themselves are inhomogeneous, as seen 
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in the detailed studies of NGC5471 by Skillman (1985), and of NGC604 by Diaz et al. (1987). That there 
are also inhomogeneities on smaller spatial scales is shown by the discrepancy between (rms) values 
derived from recombination emission and local values determined from density-sensitive line ratios. The 
dense clumps are embedded in a much lower-density medium, with typical clump volume filling factors of 
0.01 -0.1 (e.g. Kennicutt, 1984; McCall et al., 1985). The interclump material is often treated as a vacuum 
in nebular models, because it does not contribute significantly to the optical emission lines. 

Giant extragalactic Hll regions display supersonic velocities, which appear to correlate with Hj8 lumi- 
nosity. Terlevich & Melnick (1981) interpret the line- widths as virial and therefore usable for determining 
the local gravitational field; they also find a secondary dependence on metallicity. An alternative interpre- 
tation of the origin of the line-widths is that they ai^e a result of stellar winds from the exciting stars, and 
possibly also from embedded supernova remnants (e.g Dopita, 1981; Skillman, 1985). For neai^by regions, 
it is possible to actually identify the stars which may be responsible for driving the high-velocity gas. 

As mentioned above, luminous extragalactic H II regions are ionized by OB associations. For nearby 
regions, the members of the stellar cluster can be distinguished individually and HR diagrams can be con- 
structed. The nebular ionization structure and emitted spectrum will evolve as the cluster ages and the UV 
radiation field diminishes and softens. Wolf-Rayet stars are often present in extragalactic Hll regions, the 
frequency of Wolf-Rayet stars is higher for higher-metallicity regions as proved by Maeder et al. (1980). 
Wolf-Rayet stars are important in this context because they furnish metal-rich outflows which ai^e capable 
of altering the chemical composition of their gaseous environment. Giant H II regions are also known hosts 
of Type II supemovae. Winds and supernovae from the massive stars can contaminate (or enrich) the local 
gas in Hll regions in He, C, O, and other species. Evidence for such local enrichments has been sought and 
perhaps seen in some regions (Skillman, 1985; Pagel, 1986). 

1.2 Determination of chemical abundances in H n regions 

H II regions are ideal places to determine the abundance of the elements that are responsible for recombina- 
tion and fine-structure lines. The list of these elements is generally limited at the present time, although lines 
of many more elements are observable in several planetary nebulae and high spectral resolution and multi- 
wavelength studies of nearby Hll regions (e.g. Garcia-Rojas et al., 2006). The determination of element 
abundances in H II regions are given relative to the hydrogen content, which is observed by its recombina- 
tion lines. Only a few abundant elements give observable recombination lines with similar physics: helium, 
carbon, nitrogen and oxygen, whose lines are very weak to detect and may suffer problems with fluores- 
cence excitation. The abundances derived from the fine-structure lines in the visible are sensitive to both 
temperature and density, and the interpretation of the line intensities is a delicate problem. In some cases, 
the temperature of the emitting zone can be obtained and the abundance determination is safer. However, 
even in this case temperature fluctuations can yield systematic errors in the abundances. 

The abundances derived from the mid-and far-infrared fine structure lines are not sensitive to electron 
temperature and are little affected by extinction. These are considerable advantages with respect to the 
optical lines. However, there are important discrepancies between the abundances derived from infrared 
and from optical lines. These differences may originate in temperature fluctuations or in errors in atomic 
parameters, but one has to consider that the critical density for infrared lines is generally much smaller than 
for visible lines, so that the abundances derived from the infrai^ed lines ai^e underestimated if the density is 
high (this effect can be important for planetary nebulae). 
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Observed emission line intensities 

Correction for extinction using hydrogen recombination decrement 

Reddening-corrected line intensities 
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Local physical conditions (n^, TJ 

Ionic level populations and calculated line emissivities 
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CoiTection for unobserved ions 

Element abundance ratios (0/H, etc.) 

Figure 1.1: The direct method of chemical abundance determinations. 

All elements (with the obvious exception of H) exist in several ionization states in H II regions. However, 
only the abundances of those ions that emit observable lines can be determined. If such ions are minor 
species, they yield no useful information because the physical parameters of H II regions are most often 
too uncertain to allow an accurate solution of the ionization equilibrium. This is for example the case for 
O I, C II, S II or Si II. If the observed ion is a major species the situation is more favorable since we can 
calculate, more or less accurately, the abundances of the unobserved ions of the same element. However, 
uncertainties remain if a high precision is required, as is the case of helium in a cosmological context. The 
most favorable case is that of oxygen, whose major ionization states, O II and O III, are observable optically 
and for which the electron temperature Tg can be determined. For this reason, oxygen is, after helium, the 
element whose abundance is best determined, at least if the temperature is large enough (or the metallicity 
is too low) for the temperature-sensitive lines (e.g. [O III] A4363) to be measured. If this is not the case, we 
may construct tailored models of the nebular ionization and thermal structure of a H II region to estimate the 
electron temperatures and ionization correction factors for individual ions. 

However, given the difficulty of detecting the Tg-sensitive line and the assumptions made in nebular mod- 
eling, a very popular approach is to obtain the abundance of extragalactic H II regions using empirical rela- 
tions between the oxygen abundance and the intensity of the [O ll] A A 3726, 3729 and [O ill] AA4959, 5007 
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lines relative to Hj3 (Pagel, 1997) or by using the [Oil] AA7320,7330 as described by Aller (1984) and 
implemented by Kniazev et al. (2004) in SDSS Hll galaxies. This method however, is the less accurate and 
much discussion about the reliability of the different empirical calibrations is still ongoing in the Uterature 
(e.g. see Kewley & Ellison, 2008, for a thorough discussion). A full discussion regarding this topic is beyond 
the scope of this chapter, however, in ??, I include a small review on the different empirical techniques of 
abundance determinations (considering their particular advantages and pitfalls), and their implementations 
in the context of the work carried out in this dissertation. A more complete explanation of the determination 
of nebular abundances from emission lines can be found in references on the physics of gaseous nebulae 
such as Aller (1984), and Osterbrock & Ferland (2006). 

1.3 Abundance gradients in galactic disks 

It has been noticed that certain H II region emission-line ratios, such as [O III]/Hj8 , vary across the disks of 
nearby spiral galaxies. The interpretation of this variation in terms of a metallicity trend was introduced 
by Searle (1971), in a paper that laid the groundwork for the entire field of abundance gradients. It was 
soon followed up by further observational studies and a more rigorous analysis involving the construction 
of realistic nebulai^ models (Shields, 1974). From the start it was recognized that there was a need for a 
"second parameter" in addition to the 0/H ratio, to explain an observed systematic increase in 0+^/0+ with 
decreasing 0/H. Shields & Tinsley (1976) suggested that this secondary effect results from a tendency for 
the effective temperatures of the ionizing stars to be hotter for lower 0/H, and interpreted it as a metallicity- 
dependent truncation of the top end of the initial mass function (i.e. that the formation of very massive stars is 
inhibited by higher metallicity). Some form of the idea of a Z-dependent IMF is still a popular interpretation 
of the "excitation" trend (e.g. Vilchez & Pagel, 1988), but it is also the case that a similar effect can arise 
from systematic variations in the nebulai^ geometry and/or filling factor (Mathis, 1985; Dopita & Evans, 
1986). 

An extensive body of literature has been amassed on the subject of abundance gradients in galaxies. Not 
surprisingly, many works have focused on large, nearby galaxies with many observable H II regions, such 
as M 33 (Vilchez et al., 1988; Rosolowsky & Simon, 2008) and M 101 (Evans, 1986; Torres-Peimbert et al., 
1989; Kennicutt & Garnett, 1996). The gradients are usually expressed as a logarithmic fit to some 5-20 
regions per galaxy, and have a magnitude of about 5log(0/H)/R = -0.08 (± 0.03) dex/kpc. This is similar 
to the values derived for the solar-neighborhood metallicity gradient in the Milky Way galaxy. The trend of 
these gradients in the inner parts of galactic disks ai^e difficult to study, both because the H II region samples 
are often small, and more fundamentally because these are generally the most metal-rich regions, for which 
[O III] A4363 is unobservable and therefore the derived abundances are heavily model-dependent. 

The steepest abundance gradients were initially seen in late-type spiral galaxies (types Sb-Scd). Irreg- 
ulars and barred spirals tend to have weak or zero radial gradients. Early type spirals are harder to study 
because their H II regions are intrinsically fainter, but studies of M 81 (Sab) show it to have an 0/H gradient 
similar to those of M 33 and M 101 (Garnett & Shields, 1987). There is at present no convincing evidence 
that the 0/H gradient depends on morphological type among spiral galaxies. However, there is evidence for 
a good correlation between mean 0/H abundance and the overall galaxy mass or luminosity. This trend re- 
sembles the correlation of stellar metallicity with galaxy mass, and probably has its roots in the fundamental 
processes of galaxy formation and evolution. 

Along with the trend in [OlII]/Hj8, a similar radial trend was noted for the ratio [NlI]/Ho;, which de- 
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creases with increasing distance from the centers of spiral galaxies. Although part of this trend is due to the 
generally lower degree of ionization in the outer H II regions, there also must be a real variation in abun- 
dance. Unlike oxygen, for nitrogen one usually can measure the singly-ionized state only; unfortunately, 
N++ has no strong optical lines. As mentioned before, the nitrogen abundance is basically derived from 
[NlI]/[OlI]. The relative behavior of O and N is often displayed by plotting N/0 vs. 0/H. Some studies 
find that N/0 varies almost as steeply as 0/H, which has special significance in the context of chemical 
evolution models, but others claim that N/0 varies only slightly or is constant across the disks of galaxies 
such as M 101, M33, M81, and M83. There also appear to be variations in N/0 at a given 0/H from 
galaxy to galaxy (same references as above). Some of these variations may be an artifact of the analysis, 
especially since N+ contains only a small fraction of the nitrogen for the lowest-abundance, most highly ion- 
ized regions. For such regions, the ionization correction factors are very large, and the uncertainties in the 
ionization structure translate into large uncertainties in the elemental abundance of nitrogen. Nevertheless, 
there is accumulating evidence that nitrogen has a more complicated behavior than does oxygen, with N/0 
being roughly constant at low values of 0/H and increasing at higher 0/H (e.g. Pagel, 1985; Torres-Peimbert 
et al., 1989). Measurements of N/0 in metal-poor dwarf irregular galaxies are an important ingredient in 
this argument. 

Scatter in gradient determinations has been seen in various studies (e.g., in the Milky Way Afflerbach 
et al. 1997 or in M 33 Rosolowsky & Simon 2008), even after accounting for uncertainties in the stellar 
absorption and reddening corrections, an intrinsic scatter of ~ 0.1 dex has been measured in these very 
well-studied galaxies which is unexplained by the measurement uncertainties. Regardless of its source, 
gradient determinations made in the face of significant scatter coupled with a limited number of observations 
may produce widely varying results. This historical evolution of the gradient determination ranging over 
nearly an order of magnitude, should serve as a cautionary example. Only lai^ge numbers of measurements 
can overcome the uncertainties engendered by the intrinsic variance, as some observations suggest that the 
uncertainties in the gradients are systematically underreported. 
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1.3.1 The Galactic abundance gradient 

Because of interstellar extinction, one can use the same techniques as for extragalactic H II regions only 
for the part of our Galaxy outside a galactocentric distance of about 7 kpc. Studies such as those by Haw- 
ley (1978) found gradients similar to those in other spirals, 5log(0/H)/5R = -0.04 to -0.06 dex/kpc and 
5log(N/H)/5R = -0.10 dex/kpc. Determination of abundances in the inner galaxy requires the use of other 
techniques, such as measuring electron temperatures from radio recombination lines. The values of are 
found to increase systematically with increasing radius, presumably because of a decreasing abundance of 
oxygen, the primary coolant. The inferred gradient in 0/H is 5log(0/H)/5R = -0.07 dex/kpc after the classic 
paper of Shaver et al. (1983). 

The results from optical studies for the other measurable elements ai^e similar- to those for other galaxies: 
N/H varies more steeply than 0/H; S/0, Ne/0, and Ar/0 do not vary in the outer part of the Galactic disk. 
Again, the optical studies are restricted to the unobscured portion of the Milky Way galaxy, and therefore 
do not sample the inner disk where the inferred 0/H values are high. A more recent development, made 
possible by improvements in infrared detectors and the availability of space observatories. The exploration 
of the infrared spectral region as a tool for studying the galactic abundance gradient. The mid-infrared 
spectral region (5-30 /^m) contains emission lines of the major ions of Ar, S, and Ne: [Aril] 7.0 and [Arm] 
9.0 /im; [S III] 18 and [S IV] 10.5 /xm; and [Nell] 12.8 /xm. These lines have been measured in a number of 
H II regions in the inner Galaxy, and evidence for abundances elevated by factors of two or three have been 
found for the Galactic Center and for H II regions in the 5 kpc "ring" region (Pipher et al., 1984). 

However, even these mid-infrared lines suffer somewhat from extinction. In particular, the [Ar III] and 
[SiV] lines fall in the middle of the strong 10 /im silicate absorption feature, where the optical depth is 
comparable to that in the near-infrared. Another approach to studying abundances in the inner galaxy is 
to make use of the fine-structure lines of [OlII] 52, 88 jxra and [NIII] 57 /xm. By a happy coincidence, 
these hnes from the abundant and (presumably) usually co-extensive 0++ and N++ ions fall close together 
in wavelength and have fairly similar dependences on the electron density. The line emissivities are also 
essentially independent of the electron temperature. Measurements of these three lines therefore yield a 
relatively accurate value for the N/0 ratio. A survey of about a dozen galactic H II regions in these lines 
yielded strong evidence that N/0 in the Galactic Center and 5 kpc "ring" is elevated by a factor of 2 or 3 as 
compared to the solar neighborhood. There remain some unsettled questions regarding N/0 determinations 
from the far-infrared lines, including possible ionization structure effects in Hll regions ionized by very 
cool stars and a systematic discrepancy between values derived from the infrared lines and those derived 
optically from [Nll]/[Oll] (Rubin et al., 1988). 

More recent observations of IR fine-structure hnes of the [S III] 19 /im, [O III] 52 and 88 /xm, and [N III] 
57 jxra in compact H II Galactic regions have found abundance gradients of the form [S/H] = (-4.45 ± 0.04) - 
(0.063 ± 0.006) (kpc), [N/H] = (-3.58 ± 0.04) - (0.072 ± 0.006) (kpc), and [0/H] = (-2.85 +/- 0.06) - (0.064 
+/- 0.009) (kpc) (Afflerbach et al., 1997). These abundances are consistent with production of sulphur, 
nitrogen, and oxygen by primary nucleosynthesis. Comparison with abundances in other galaxies implies a 
Hubble type between Sab and Sb for our Galaxy and an unbaiTcd or mixed galactic structure (Vila-Costas 
& Edmunds, 1992). 
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1.4 Comparison with Chemical Evolution Models 

The recognition of significant variations in the gas composition within and among galaxies, along with 
parallel results on the stellar populations, inspired the development of chemical evolution models which 
attempt to explain these patterns. The so-called "simple model" postulates a closed system of gas and stars, 
which self-enriches in metals as generations of stars age, die, and seed the ambient gas in the heavy elements 
(Seai^le & Sargent, 1972). This model also makes the approximations that the stellar lifetimes and timescale 
for complete mixing of nucleosynthetic products are negligible in comparison to the timescale on which the 
metallicity evolves ("instantaneous recycling"). The simple model makes a specific prediction regarding the 
metallicity and system properties: 

Z = yln{Mtota!/Mgas), (1.1) 

in this equation Z is the metal abundance, y is the fraction of the stellar mass converted to heavy elements 
(the yield), and Mfotai = Mgas + M, 

tars- 

Although this model is most appropriate for the low-mass galaxies, it can also be applied to large disk 
galaxies if concentric radii are treated as independent zones. However, it does not explain the observed gra- 
dients, so modifications such as radial flows, matter exchange with an outside reservoir (infall and outflow), 
or a variable stellar initial mass function, have been proposed as modifications to the model (Matteucci & 
Francois, 1989; Dopita, 1990). 

The relative abundances of nitrogen and oxygen are of particular interest, since they are synthesized 
in different astrophysical sites. Oxygen is synthesized in massive stars and distributed into the interstellar 
medium by Type II supemovae, while the origin of nitrogen is more problematic. A distinction is frequently 
made between "primary" nucleosynthetic products, which can be synthesized directly from H and He in 
Population III stars, and "secondary" products, which require a "seed" heavy nucleus to be initially present 
in the star where its synthesis occurs. By this definition, oxygen is a primary species. Nitrogen is secondary 
when made as a by-product of CNO-cycle hydrogen burning. According to the simple closed-box model, 
the abundance of a secondary species is quadratic, so that if N is secondary and O primary, then (N/H) 
oc (0/H)^, or (N/0) (0/H). The N/0 ratio does appear to approach this behavior, for Hll regions with 
moderately high 0/H values in M 101 (Torres-Peimbert et al., 1989) and in the Milky Way. However, below 
a certain values of 0/H, it appears that N/0 is constant; these low-metallicity Hll regions occur mostly 
in low-mass galaxies. Thus, it is becoming clear that nitrogen is not purely a secondary nucleosynthetic 
product. Indeed, N may be produced within intermediate-mass stars by an effectively primary process, if C 
synthesized within the star by the triple-alpha reaction is later subjected to the CN cycle. Nitrogen made by 
this process would be primary, but there might be a time-delay in building up its abundance relative to the 
nuclear products of supemovae, because of the longer lifetimes of the source stars. 

The other elements measured in extragalactic Hll regions, S, Ne, and Ar, are not likely to be dominated 
by secondary processes. They might still, however, vary differently than oxygen, if they were produced in 
stars of different mass ranges and the IFM varied or the timescales for enrichment differed substantially. 
There are known variations in the abundance ratios of certain elements. For example the fact that the iron- 
group is deficient relative to oxygen in Population II stars is thought to reflect an origin for the former 
chiefly in Type I supemovae, which originate in long-lived progenitors, as opposed to synthesis of oxygen 
in massive stars and Type II supemovae. 
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In the context of chemical evolution models, Garnett (2002) studied the metallicity-luminosity and 
metallicity-rotation speed correlations for spiral and in^egulai^ galaxies for a sample of spiral and irregu- 
lar galaxies having well-measured abundance profiles, distances, and rotation speeds. He finds that the 
O/H-Vro, relation shows a change in slope at a rotation speed of about 125 kms^'. At faster there 
appears to be no relation between average metallicity and rotation speed. At lower "Vrot, the metallicity 
con^elates with rotation speed. This change in behavior could be the result of increasing loss of metals from 
the smaller galaxies in supernova-driven winds. The idea was tested by looking at the variation in effective 
yield, derived from observed abundances and gas fractions assuming closed box chemical evolution. The 
effective yields derived for spiral and in^egular galaxies increase by a factor of 10-20 from V,-o, ~ 5 to 300 
km s~^ asymptotically increasing to approximately constant ye/f for ~ 150 km s^^. The trend suggests 
that galaxies with ~ 100-150 kms^^ may lose a large fraction of their supernova ejecta, while galaxies 
above this value tend to retain metals. The determination of effective yields as function of galactic radius 
and its interpretation stands as one of the main studies in order to discriminate among different physical 
effects which may affect the chemical evolution of a galaxy. 

1.5 Goals of this dissertation 

As described in this chapter, the study of chemical abundances has undergone a remarkable development 
in the last decades thanks mostly to important observational efforts that have focused on the derivation of 
physical and chemical properties of emission line Hll regions in galaxies by spectroscopic techniques. The 
main motivation common to all of these observations is to use the chemical information as one of the means 
at our disposal to understand the physical processes at play in the formation and evolution of galaxies in the 
universe. 

Hitherto, most spectroscopic studies in nearby objects have been limited by the number of objects sam- 
pled, the number of H II regions observed and the coverage of these regions within the galaxy surface. In 
order to increase significantly the number of H II regions sampled in any given galaxy we require the com- 
bination of high quality multi-wavelength data and wide field spectroscopy. The advent of multi-object 
and integral field spectrometers with lai^ge field of view now offer us the opportunity to undertake a new 
generation of observations, based on samples of scores to hundreds of H II regions and full 2-dimensional 
(2D) coverage. These sort of data would enable to test, confirm and extent the previous body of results from 
small sample studies, while at the same time open up a new frontier of studying the 2D metallicity structure 
of disks and the intrinsic dispersion in metallicity, or to test and strengthen the diagnostic methods that ai"e 
used to measure the Hll region abundances in galaxies, among other issues. 

The scientific core of this dissertation is based on an observational project conceived to tackle the prob- 
lem of the 2D spectroscopic coverage of the whole galaxy surface. New techniques in imaging spectroscopy 
(or integral field spectroscopy, IFS) provide a powerful tool for studying the small and intermediate scale- 
size variation in line emission and stellar continuum in nearby well-resolved galaxies. We designed a project 
to take advantage of these new observational techniques in order to assemble a unique spectroscopic sam- 
ple from which we could study, with unprecedented detail, the star formation and gas chemistry across 
the surface of a galaxy. The observations consist of Integral Field Unit (IFU) 2D spectroscopic mosaics 
of a representative sample of neai^by galaxies (D < 100 Mpc) with a projected angular size of less than 
10 arcmin. The mosaics were constructed using the unique instrumental capabilities of the Postdam Multi 
Aperture Spectrograph, PMAS (Roth et al, 2005) in the PPAK mode (Verheijen et al, 2004; Kelz & Roth, 
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2006) at the German-Hispanic Astronomical Centre at Calar Alto (CAHA), Spain. The PMAS fibre PAcK 
(PPAK) is one of the world's widest integral field unit with a field-of-view (FOV) of 74 x 65 arcseconds that 
provides a semi-contiguous regular sampling of extended astronomical objects. This project represents the 
first attempt to obtain 2D spectra of the whole surface of a galaxy in the nearby universe. The spectroscopic 
mosaicing comprises more than 50 000 spectra in the optical wavelength range. 

This observational project was devised as a scientific international consortium, the members ai^e world- 
leading experts in their respective fields, including star formation and chemical abundances of galaxies, 
active galactic nuclei, multiwavelength observations of emission line regions and 2D spectroscopy. The 
project was entitled: the PPAK IFS Nearby Galaxies Survey, or PINGS. The P.I. of this project is Prof. 
Robert C. Kennicutt Jr. at the Institute of Astronomy, University of Cambridge. The collaborators of the 
consortium are: Dr. Angeles Diaz at the Universidad Autonoma de Madrid, Spain; Dr. Anna PasquaU at 
the Max-Plank Institut fiir Astronomie in Heidelberg, Germany; Dr. Sebastian S. Sanchez at CAHA, Spain; 
Benjamin Johnson and Caina Hao at the University of Cambridge, UK. 

The primary scientific objectives of this dissertation are to use the PINGS observations to obtain pixel- 
resolved emission-line maps across the disks of the galaxies to study the 2D abundance distribution and on 
characterising the relations between these abundance properties and the physical properties of the parent 
galaxies. By targeting virtually every Hll region in the galaxies, as a consequence of the nearly complete 
spatial coverage of the IFUs, we are able to test for the first time the systematic dependences of the strong- 
line abundances on the size, luminosity, surface brightness, and other properties of the Hll regions. In 
that respect, the PINGS observations and the subsequent analysis represent a leading leap in the study 
of the chemical abundances and the global properties of galaxies, information which is most relevant for 
interpreting observations at all redshift sources accessible with the current technology. 

1.5.1 Structure of the dissertation 

The structure of this thesis is as follows: In § 2, 1 discuss the importance of Integral Field Spectroscopy (IFS) 
in astrophysics, including an explanation of the technique with their advantages and pitfalls, a description of 
the available instrumentation in the world by the time this project was envisaged, and the selection criteria 
of the telescope-instrument chosen for this project. I also include a brief review to previous works that have 
attempted to obtain 3D chemical abundance information in galaxies. In § 3, 1 present the sample of galaxies 
and the selection criteria followed according to the scientific objectives established for this dissertation. This 
chapter also includes a full description of the logistics and observations, explaining the telescope set-up and 
the particular observing technique adopted for this project. In § 4, 1 explain the reduction process of the IFS 
raw data, the additional corrections implemented in this project and the improvements with respect to pre- 
vious pipelines, particularly regarding the flux calibration and the sky-subtraction. The possible sources of 
eiTors and uncertainties are addressed, together with an explanation of the techniques implemented to min- 
imise them. In § 5, 1 present the integrated spectra of the PINGS sample, obtained by co-adding the spectra 
from their corresponding mosaics. Comparisons with previously published data are included. An analysis 
of the ionized gas component is performed, together with the techniques and methodologies implemented in 
order to derive the physical parameters of the integrated gas-phase of the galaxies sample. In § 6, 1 present 
a complete 2D spectroscopic study for a selected number of galaxies from the sample. Selected H II regions 
previously observed ai^e compai^ed with spectra extracted from the PINGS sample. A set of emission line 
maps calculated from each galaxy is presented, including a quantitative description of the 2D distribution of 
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the physical properties inferred from them. Then, a detailed, spatially -resolved spectroscopic analysis of the 
selected galaxies is performed, based on different spectral samples extracted from the full IFS mosaics of 
the galaxies. Several diagnostic diagrams and the state-of-the-art abundance diagnostic techniques are used 
to obtain the 2D distribution of the physical properties and chemical abundances of the selected sample. 
Finally, in § 7 I present the general conclusions of this dissertation, including some planned paths of future 
investigation. 
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Thesis content 



The whole thesis is not included in astro-ph due to file size limitations. 

The full contents can be found at: 
http : / / www . dspace . cam . ac . uk/handle/ 1810/224843 
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